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ABSTRACT 

In the two previous papers of this series, we have discussed the importance 
of the /^-gradients due to hehum setthng on rotation-induced mixing, first in 
an approximate analytical way, second in a 2D numerical simulation. We have 
found that, for slowly rotating low mass stars, a process of "creeping paralysis" in 
which the circulation and the diffusion are nearly frozen may take place below the 
convective zone. Here we apply this theory to the case of lithium and beryllium 
in galactic clusters and specially the Hyades. We take into account the rotational 
braking with rotation velocities adjusted to the present observations. We find 
that two different cells of meridional circulation appear on the hot side of the 
"lithium dip" and that the "creeping paralysis" process occurs, not directly below 
the convective zone, but deeper inside the radiative zone, at the top of the second 
cell. As a consequence, the two cells are disconnected, which may be the basic 
reason for the lithium increase with effective temperature on this side of the dip. 
On the cool side, there is just one cell of circulation and the paralysis has not 
yet set down at the age of the Hyades; the same modelisation accounts nicely for 
the beryllium observations as well as for the lithium ones. 

Subject headings: stars: abundances; stars: rotation-induced mixing; diffusion 

1. Introduction 

This paper is the third of a series on the coupling between helium settling and meridional 
circulation in stellar radiative zones. The process of meridional circulation and the way it 
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can slow down diffusion have been studied extensively in the past by many authors (see 
Vauclair & Theado (2002), hereafter referred as paper I, for references). From the beginning 
of the computations of meridional circulation and rotation-induced mixing, the importance 
of the feed-back effect due to nuclearly-induced ^-gradients was well recognised (see Mestel 
& Moss 1986), however the feed-back effect due to the diffusion-induced //-gradients on the 
meridional circulation was not included until recently. Vauclair (1999) showed how, in slowly 
rotating low mass stars, the resulting terms in the computations of the meridional velocity 
could become of the same order of magnitude as the other terms. Theado & Vauclair (2001) 
computed these terms for the case of Pop II stars with the assumption of negligible differential 
rotation. They claimed that the induced diffusion-circulation coupling could be the reason 
for the very small dispersion of the lithium observed in these stars. More recently (paper I) 
the effect of differential rotation was discussed and we showed how it could be introduced 
in a general context. Applications to different stellar populations are underway. Here we 
discuss the case of Pop II main sequence stars. 

While the /x-terms have definitely to be taken into account in the computations, the 
details of the coupling process which occurs when these terms becomes of the same order 
of magnitude as the classical terms remain difficult to handle. In paper I, we introduced 
an analytical approach of the processes and discussed an approximate solution in a quasi- 
stationary case. We showed that the couphng between hehum settling and rotation-induced 
mixing can slow down both the mixing and the settling when the //-terms are of the same 
order of magnitude as the other terms. In paper II (Theado & Vauclair 2002), we presented 
the results of a 2D numerical simulation of the considered processes which helps visualizing 
the situation. In the example which has been chosen (slowly rotating low mass stars), the 
/i-terms become of the same order of magnitude as the classical terms, in a short time 
scale compared to the main-sequence lifetime, thereby creating a frozen region where the 
circulation does not proceed anymore. However when this occurs the equilibrium between 
the opposite meridional currents is permanently destabilized by the helium settling below 
the convective zone. As a consequence, a new motion develops, which mixes up the zone 
poUuated by diffusion. Below this region the circulation remains nearly frozen. This process 
acts as a restoring system for the element abundance variations and increases in a significant 
way their settling time scales. We then suggested that diffusion and mixing react in such a 
way as to keep both the horizontal and vertical /t-gradients constant in the frozen region, 
while they proceed freely below. 

We now apply this prescription to light elements in population I stars. In this paper, we 
treat the case of the main-sequence F and G stars. In forthcoming papers we will extend our 
study to the case of sub-giant F stars and we will discuss the effect of differential rotation 
in the case of Pop II stars. We emphasize that up to now these studies have been applied 
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to low mass stars with deep convective zones, in which the Gratton-Opik term (see eq. 3) 
is neghgible. This is not the case for F stars with effective temperatures larger than about 
6700K (hot side of the lithium dip in galactic clusters). Wc will show that in this case 
two different loops of meridional circulation may arise and that the freezing effect on the 
circulation does not take place in the first loop, below the convective zone. It may however 
have some influence at the boundary between the two zones, as discussed below. 

In section 2, we review the observations of lithium abundances and rotation velocities 
in the Hyades. 

A theoretical discussion about previous studies is given in section 3. 

In section 4 we recall some important expressions of the meridional circulation and the 
rotation induced-mixing. 

In section 5 we present the results of the evolutionary computations. We show the 
influence of the feed-back effect due to the diffusion-induced yU-gradients on the meridional 
circulation. We discuss the important difference in the implied physical processes which 
occur on the hot and cool sides of the lithium dip. 

In section 6, we compare the results of the computations with the lithium and beryllium 
observations in galactic clusters. We give our conclusions and discuss future work. 

2. Observational constraints 

2.1. The lithium dip 

Lithium abundance determinations in population I stars show that lithium is strongly 
depleted in the main-sequence stars located in a narrow range of 300K in effective temper- 
ature around 6700K. This feature called the "lithium dip" had been observed for the first 
time in the Hyades by Wallerstein et al. (1965) and much later by Boesgaard & Tripicco. 
(1986). Since these first observations many abundance determinations have confirmed the 
presence of the dip in all galactic clusters older than 10^ yrs (NGC 6475 (0.2 10^ yrs), Ursa 
Mayor (0.5 10^ yrs), Pracscpc (0.8 10^ yrs). Coma Berenice (0.8 10^ yrs), NGC 752 (2.2 10^ 
yrs)) as well as in field stars. In younger clusters (Alpha Persei (0.05 10^ yrs), Pleiades (0.08 
10^ yrs)), the phenomenon is not observed, although a small decrease in lithium abundances 
is observed in the Pleiades around 6700K. In clusters older than NGC 752 (M67 (5 lO'^ yrs), 
NGC 188 (6.5 10^ yrs)), the stars of the lithium dip already evolved off the main-sequence. 

On the red side of the dip, the lithium abundances increase for smaller effective tem- 



-4- 



peratures until they reach a sort of plateau, which extends from 6400 to 6000K. In cooler 
stars (G type stars) the lithium abundances decrease with effective temperature. 

On the blue side of the dip, the lithium abundances increase with the effective temper- 
ature. Few abundance determinations are available for main-sequence cluster stars on this 
side of the dip because most of these stars are fast rotators whereas lithium can be observed 
in slow rotators only. For Teff greater than 6900 K, the lithium abundances seem to remain 
constant, close to the galactic value (Log N(Li) ~ 3.31 with Log N(H)=12); only few Am 
stars display anomalous abundances. 

Lithium observations in the subgiant stars which originate from the hot side of the 
dip bring some interesting constraints on hydrodynamical processes which occur during the 
main-sequence phase. The convective envelopes of the subgiants and giants have deepened 
and dredged-up material from the interior of the stars and thus provide a view of the internal 
composition. Many observations show that subgiant stars originating from this side of the dip 
present large range of lithium abundances which cannot be explained by standard dilution 
alone (Alschuler 1975, Brown et al. 1989, Balachandran 1990, Dias et al. 1999). These 
observations suggest that a lithium depletion process takes place in the radiative interior 
during the main-sequence phase even if its effects are not visible at the surface at the age of 
the Hyades (Vauclair 1991, Charbonnel & Vauclair 1992). 

2.2. Rotational velocities in the Hyades 

Boesgaard (1987) noticed that the hthium dip coincides in the Hyades with a sharp 
drop in rotational velocities. The determination of rotational velocities in the Hyades (Gaige 
1993), Praesepe and Coma Berenice clusters shows that stars cooler than 6000K are slow 
rotators (excepted a few M stars). The increase of the rotational velocity coincides with the 
lithium decrease on the red side of the dip. On the blue side, the influence of the velocity is 
more difficult to determine because of the observational bias previously mentioned : in hottest 
stars {Teff > 7000 K), the lithium line is very weak and is detectable only in stars whose 
rotational velocities are small (Vrot < 50 km.s~^). When the effective temperature decreases, 
lithium can be observed in stars with higher rotational velocities. This observational bias 
induces a second bias since slow rotator with Teff > 7000K are chemically pecuhar stars. 
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3. Previous theoretical studies 

Many explanations have already been proposed to account for the observed lithium 
abundances in the Hyades. Most of them rely on the fact that the lithium dip appears in the 
special range of effective temperature where there is a transition between stars with deep 
convective zones (on the cool side) and shallow ones (on the hot side) . We will here briefly 
summarize them. 

The earliest and simplest explanation for the lithium dip was the microscopic diffusion 
model of Michaud (1986) (see also Proffitt & Michaud (1991), Richer, Michaud and Proffitt 
(1992) and Richer & Michaud (1993) for more sophisticated treatments of the diffusion). 
In the temperature range of the dip the convective depth rapidly decreases with increasing 
effective temperature which strongly affects the microscopic diffusion. They predicted that 
lithium would settle gravitationally out of the convective zone in stars with Teff < 6900K. 
At hotter temperatures, upward radiative acceleration was predicted to dominate over grav- 
itation, leading to surface lithium enrichment which are not observed. To explain the lack 
of lithium overabundances, Michaud et al. invoked mass loss. The rise of the lithium abun- 
dance on the red side of the dip was explained by the increase of the convective zone depth : 
in stars cooler than 6400K, it leads to longer diffusion time scales and thus to slower surface 
lithium depletion. They finally succeeded in reproducing the lithium dip by adjusting two 
free parameters : the mixing length and the mass loss rate. 

However several observational facts argue against the pure microscopic diffusion hy- 
pothesis. First the predicted dip is narrower than observed (Balachandran 1995). Second 
the carbon, oxygen and boron under-abundances expected in the case of pure diffusion are 
not observed in the Hyades (Boesgaard 1989, Priel & Boesgaard 1990, Garcia Lopez et al. 
1993). Finally, according to the models, gravitationally settled Li is supported in a region 
just below the convective zone and is dredged up to the surface as soon as the star evolves 
off the main-sequence. Observations on M67 (Balachandran 1995) show that this lithium 
dredge-up is not observed at the turn-off or on the sub-giant branch. 

Based on the rough coincidence between the instability strip and the position of the 
Li dip, Schramm et al. (1990) suggested that pulsation-driven mass loss could account 
for the lithium dip. This hypothesis is in severe disagreement with several observations 
published since their study. It would mean that no lithium depletion should be perceived 
in the mass-losing model until the convective zone reaches the top of the lithium nuclear 
destruction region, at which time the depletion would be extremely rapid : the observations 
which indicate that the dip forms gradually contradict this result. Moreover Delyannis & 
Pinsonneault (1993) found that beryllium was also depleted in these stars, in contradiction 
with this theory : according to the mass loss model the hthium should be totally destroyed 
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before Be depletion appears. 

Garcia Lopez et al. (1991) argued that internal waves generated by the convective 
envelope of a star may be effective in producing a weak mixing in the radiative interior. 
They succeeded in reproducing the lithium dip but only by increasing the intensity of the 
waves by a factor 15 above the estimates obtained with straightforward mixing length theory. 
Moreover their models fail to reproduce the large abundance dispersion on the red side of 
the dip. 

Many theoretical works have also tested the influence of the rotation-induced mixing on 
the dip stars. Charbonnel et al. (1992) (hereafter CVZ92) studied the lithium depletion by 
meridional circulation-induced mixing. They used a simplified expression of the meridional 
velocity which did not take into account the /i-gradicnts neither the differential rotation. In 
this case, the meridional velocity vanishes at the depth where p = VL^/{2'kG). Below this 
layer the circulation goes in the classical sense, while it is inverted above : the circulation 
is divided in two separate cells. CVZ92 pointed out that the separation of the circulation 
in two cells occurs precisely at the effective temperature of the bottom of the lithium dip: 
the second cell appears on the hot side while the boundary is inside the convective zone 
on the cool side. This result helped them give a consistent modelisation of the dip. Under 
the effects of the circulation-induced mixing they showed that the theoretical predictions 
and the observations agreed on the cool side of the dip. For the hot side of the dip they 
invoked the "quiet zone" which separates the two cells of meridional circulation. However it 
was later pointed out (Michaud, private discussion) that this zone was too thin to prevent 
microscopic diffusion so that lithium depletion subsisted anyway. In the present paper, we 
will revisit this theory in the new framework of rotation-induced mixing including /x-gradients 
and differential rotation. 

Talon & Charbonnel (1998) used the lithium observations in the Hyades to assess the 
role of meridional circulation and shear turbulence in the transport of angular momentum. 
Rotation-induced mixing was included as described by Zahn (1992). They computed simul- 
taneously the transport of chemicals and the transport of angular momentum due to the 
wind-driven meridional circulation. In this framework, they reproduced the hot side of the 
lithium dip but not the cool side. As this prescription also fails to reproduce the fiat rotation 
profile in the sun, they deduced that the angular momentum transport is different on both 
sides of the dip : it would be due to rotation-induced mixing on the blue side, and to another 
still unknown process on the red side. 

When the /i-gradients are not taken into account, the two loops of circulation do not 
appear on the hot side of the dip. When they are introduced, the two loops reappear 
(Palacios et al 2002) : this will be discussed again below. 
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4. Diffusion and mixing in the presence of //-gradients 

4.1. Rotation-induced mixing 

Here we recall some important expressions of the rotation-induced mixing. Following 
papers I and II, the vertical velocity amplitude of the meridional circulation can be written: 
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• the term related to the time variations of the differential rotation 
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The last term in the expression of is derived from the last term of equation (4.20) 
in Maeder & Zahn (1998) (hereafter MZ98); in the following of their study, MZ98 chose to 
neglect it; here this term is kept in the equations because it can become important in stars 
more massive than 1.3Mq. Note also that in MZ98, T which appears here in the expression 
of Ei^, was forgotten in equation (4.38), as well as in equation (4.36). 

The fourth term in equation (2) has a special behavior. It appears in the case of 
large horizontal turbulence (see MZ98), due to the fact that it modifies the radiative energy 
transport. This term can be written : 



Eh — -^—^CpTDhC 



(6) 



or 
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where K is the thermal diffusivity and D/j the horizontal diffusion coefficient. 
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In the assumption of a rotation-induced transport of angular momentum, the horizontal 
diffusion coefficient is computed accordingly, although an unknown coefficient of order one 
still remains (Talon and Charbonnel 1998). In all cases, Dh depends on Ur so that this term 
has to be treated in a different way than the three other ones. 

In paper I, we have shown that C can only vary between (no differential rotation) and 
A (maximum differential rotation). If = 0, the corresponding terms disappear, while for 
the other extreme, where C = A, we can write : 



Ec 



2'kGp 

Pra 

p \ 3 dr 



9_ 

9 

TT A 

Ht^ XT A 

ar 



— ^A + -A 
r 3 



27rG'p 



A 



Er 



(8) 

(9) 
(10) 



In the stationary case the time derivative disappears (£'^=0) and, according to Chaboyer 
& Zahn (1992) : 

A=-j^f/.^ (11) 
6Dh or 
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In the present paper, the horizontal diffusion coefficient is approximated following MZ98: 
Dh = Chi^\Ur\, where Ch is an unknown coefficient. The absolute value {\Ur\) does not 
appear in MZ98 but it is here necessary since the meridional velocity amplitude can become 
negative as we will see later. Taking into account the parametrisation of D^, we obtain : 

^" 6Cn\Ur\ dr 

Ur and A have logically the same sign. A positive value of Ur (i.e. an upward flow near the 
rotation axis and a downward flow near the equator), induce molecular weights greater at 
the rotation axis than at the equator so A is positive. Inversely a negative value of Ur, as it 
will be the case below the convective zone in the stars more massive than LSM©, leads to 
molecular weights smaller at the rotation axis than at the equator and so to negative value 
of A. 

Using equation 6 (with ^ = A) and taking into account the expression of A (eq. 11) 
leads to : 

E, = -^C,fU.^-^ (13) 

Eh is then proportional to Ur while the dependence on the horizontal diffusion coefficient 
Dfi disappears. In other words, as soon as the differential rotation has the maximum pos- 
sible value (C = A), the results do not depend anymore on the kind of angular momentum 
transport which occurs in the star. 

The Eh term must then be introduced in the computations in a different way, on the 
left side of equation 1. We can write : 

,,rC.(V.^- V + V.) ^^ ^ ^ _ ^^^gM^^ ^^^^ 

and then : 

P d In Z2 d In 11 

Note that — — - — = Hp — - — = — Vu. The final expression for Ur is : 
pg or or 

The rotational mixing is then introduced in the computations following Zahn (1992) 
and MZ98 : 
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D„ represents the vertical part of the shear-induced anisotropic turbulence due to the 
transport of angular momentum. Its efficiency depends on the Richardson criterion (Maeder 
1995, Talon and Zahn 1997). Here we chose to parametrize this coefficient in terms oi Ur ■ r 
as for the other diffusion coefficients, so that we write : 



or simply : 
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where aturb = C„ + is an unknown parameter. 
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4.2. Microscopic diffusion 

The microscopic diffusion of helium and other elements is computed as described in 
paper 11. At the considered effective temperatures the radiative force on lithium may be 
important and we have to add it in the lithium diffusion computations. Here we use the 
approximate formula developped by Michaud et al. (1976) for the radiative accelerations on 
unsaturated elements : 

«,^1.7 10»^ (20) 

where Te4 = {Tgff.lO~'^) K, A is the atomic mass of lithium and R the stellar radius. 



5. Computational results 

5.1. Chciracteristics of the models 

We have made computations for 10 stellar masses : 1.0 Mq, 1.1 Mq, 1.2 Mq, 1.25 
M0, 1.3 Mq, 1.35 Mq, 1.4 Mq, 1.45 Mq, 1.5 Mq and 1.55 Mq. Tables 1 and 2 give the 
main characteristics of these models at different ages, together with the rotation velocities 
obtained as discussed below. 

We use the statistical study of rotation velocities in the Hyades performed by Gaigc 
(1993) in order to estimate the spin down associated to the stars of different masses : we 
take an initial velocity of 100 krn.s"^ that corresponds to the mean velocity of hot stars 
(Te// >7000 K). The resulting velocity at the age of the Hyades corresponds to the average 
values for stars of a given effective temperature and these velocities are obtained by using 
the Skumanich law to reproduce the spin down of the stars. 
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Table 1: Main characteristics of the models with M < 1.25Mq. 
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Table 2: Main characteristics of the models with M > 1.3M( 
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The models are computed using the solar mixing length parameter and the metallicity 
of the Hyades : [Fe/H]=0.12 (Cayrel de Strobel et al. 1997). In each model, the parameters 
Ch and aturb (equation 19) are adjusted to obtain the right lithium depletion at the age of 
the Hyades. It is interesting to note that, while aturh is nearly constant (it slowly increases 
from 0.15 for 1.5 stars to 0.3 for 1.1 stars), the Ch values needed to account for 
the observations increase from 5 for 1.5 stars up to 1300 for 1.1 stars. This is an 
indication that the horizontal diffusion coefficient must be larger than the values obtained 
with Zahn 1992 prescription. The case of the 1.0 Mq star is still different : the values of 
the coefficients needed to explain the lithium depletion are considerably larger, which means 
that another process must occur, related to the fact that the bottom of the convective zone 
is then very close to the lithium destruction layer (e.g. a tachocline mixing). This will be 
specially discussed in a forthcoming paper (Theado, Richard and Vauclair 2003, hereafter 
TRV03). 

We will now present separately the results obtained for stars with masses lower or greater 
than 1.3 M©. In low mass stars, the Eq, term keeps the same sign in the whole radiative 
interior, inducing therefore only one circulation loop. For masses greater than 1.3 Mq, the 
Gratton-Opik term included in is negative just below the convective zone and positive 
in the deep interior. This leads to a change in the sign of Ur and to two different circulation 
loops. 



5.2. The case of cool stars (Tg// <6500K, M<1.25M0) 

First we present the results obtained for the 4 coolest stars : 1.0 Mq, 1.1 Mq, 1.2 Mq 
and 1.25 Mq. 

Figure 1 to 4 display the \Eq\ and \E^j\ profiles below the convective zone at differ- 
ent evolutionary stages. Each star arrives on the main sequence with nearly homogeneous 
composition. At that time \E^\ is smaller than l^'j^l everywhere in the star. Meridional 
circulation and element diffusion can take place in the stellar radiative region. As a conse- 
quence of helium settling the /^-currents increase below the convective zone which reduces 
the rotational mixing. 

The construction of /i-gradients and the increase of /^-currents during the stellar life 
depends on the competition between several processes. 

i) For increasing effective temperature, the convective depth decreases : the microscopic 
diffusion time scale consequently decreases which tends to speed up the construction of /i- 
gradients below the convective zone (since A is strongly related to the vertical /^-gradients) . 
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Fig. 1. — Evolution of /i-currents with radius inside a M = l.OM© star. The graphs show 
the variation with depth of both and E^^. 




Fig. 2. — Same figures as Fig. 1 for a M = I.IMq star. 
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Fig. 3. — Same figures as Fig. 1 for a M = 1.2Mq star. 




Fig. 4.— 



Same figures as Fig. 1 for a M 



= 1.25M<D star. 
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ii) For increasing effective temperature, the rotational velocity increases : the Q-currents 

are therefore more efficient for increasing masses. This can lead to different effects : a) At the 
beginning of the stellar life, the advection by circulation which only depends on fi-currents 
is more efficient, the //-gradients arc then rapidly built which leads to important ^-currents. 
These currents reduce the turbulence and the mixing, b) The circulation-induced turbulence 
is very efficient at the beginning of the stellar life. It tends to homogenize the horizontal 
layer and to delay the construction of /^-gradients. 

In these models, the "creeping paralysis" which takes place when \E^\ = \Eq\ does not 
occur yet at the age of the Hyadcs, although it is nearly reached below the convcctivc zone. 
We will discuss in TRV03 the consequences which appear later in stellar evolution, for the 
special case of the Sun. 

Figure 5 presents the lithium and beryllium abundance variations with time in the 
outer layers of these four stars. They are the result of the combined effects of diffusion and 
rotation-induced mixing. We have evoked before the influence of the convective width on 
microscopic diffusion but the depth of the convective zone plays another important role in 
the Li and Be depletion rate : the distance between the bottom of the convective zone and 
the lithium and beryllium nuclear destruction layers decreases with decreasing temperature 
so these two elements are more easily depleted in the cooler stars than in the hotter. 

Figure 6 and 7 displays as an example the lithium profiles inside the 4 stars at different 
evolutionary stages. 

5.3. The case of hotter stars (Tg// >6500K, M>1.3M0) 

We now present the result obtained for the stars hotter than 6500 K : 1.3 Mq, 1.35 M©, 
1.4 Mq, 1.45 Mq, 1.5 Mq and 1.5 Mq. In these stars the Gratton-Opik term included in Eq, 
is smaller than one in the deep radiative interior and larger than one below the convective 
zone. This leads to a change in the sign of Ur- 

Figure 8 to 11 show the \Ea\ and \Ef^\ proffies in four of the considered stars below the 
convective zone at different evolutionary stages. The drop in the \Eci\ profiles corresponds 
to the place where this term changes sign. Two circulation loops then appears : the deeper 
one goes from the rotation axis to the equator, the upper one goes in the opposite direction 
from the equator to the rotation axis (see figure 12). Then the drop in the \Eq\ proffies also 
indicates the bottom of the upper circulation loop. 

As said before, in the upper circulation loop, the circulation goes from the equator to 
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the pole. Then it brings //-enriched matter to the equator and //-depleted matter to the pole. 
This leads to negative horizontal //-gradients (A < 0) and then to /t-currents which go from 
the pole to the equator (i.e : Eq < and > 0). 

As for cooler stars, each hotter star arrives on the main sequence with nearly homoge- 
neous composition. At that time is smaller than \Eq\ everywhere in the star. Under 
the combined effects of meridional circulation and element diffusion the /i-currcnts increase 
rapidly. Figure 8 to 11 show that they become of the first order in the transition layer 
between the two circulation loops. In the 1.3 Mq and 1.35 Mq stars, the region where 
\Eq\ ~ \Efj_\ is very thin while it becomes significant in the 1.40 Mq and 1.5 Mq stars : 
about 1% of the radius in the 1.4 Mq and more than 3% of the radius in the 1.5 Mq at the 
age of the Hyades. 

This leads to the settling in the star of the quasi-equilibrium stage described in paper I 
and II. The frozen region expands slowly inside the star. Between the base of the convective 
zone and the top of the frozen region, stellar matter is mixed as well as below the bottom 
of the frozen region in the radiative interior. The stellar matter in the upper and the lower 
circulation loops are then disconnected. 

The construction of /i-gradients occurs differently here than in cooler stars. A rotational 
velocity increase leads to more vigourous Jl-currents and then to a more vigourous mixing. 
In the 1.3 and 1.35 Mq stars, the mixing prevents the /(-currents of becoming of the same 
order of magnitude as Q-currents in an extended region. In the 1.4 Mq and 1.5 Mq stars, 
the diffusion becomes more efficient than the mixing and forces the /(-currents to become 
more important. 

The lithium abundance variations with time have been computed using the following 
scenario : microscopic diffusion damped by rotation-induced mixing occurs normally until 
the relative difference between \Ef^\ and \Eq\ becomes smaller than 10%. Then, in the 
concerned region, mixing is strongly reduced while settling proceeds but with the constraint 
of keeping constant /i-gradients. Below and above the frozen region, diffusion and mixing 
are assumed to proceed freely. 

Figure 13 presents the lithium and beryllium abundance variations with time for the 
four considered stellar masses. In the 1.3 Mq and 1.35 Mq stars, the frozen region is not 
thick enough to prevent the combined effects of the microscopic diffusion and the mixing 
: lithium is depleted. In the two other stars, the frozen region is larger and significantly 
reduces the lithium depletion. The more efficient the diffusion, the more rapid the /i-currents 
increase. It explains why this process settles down earlier in the 1.5 Mq star than in the 1.4 
Mq star and why the lithium depletion is less important in the more massive star. Figure 14 
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and 15 shows the hthium and beryUium profiles in the star at different evolutionary stages. 

5.4. Discussion 

The result obtained in this paper are compared to the observational points given by 
Boesgaard & King (2002), Thorburn et al. (1993), Boesgaard & Budge (1988), Boesgaard 
& Tripicco (1986). The original lithium and beryllium abundances are the meteoritic ones 
: A(Li)=3.31 dex and A(Be)=1.42 dex, both values are from Anders & Grevesse (1989). 
Figure 16 shows the obtained results. The open symbols represent the observational points. 
The solid ones represent our models. 

All the explanations which have been proposed up to now to account for the "lithium 
dip" in the Hyades have arbitrary parameters which are adjusted at some point. For this 

reason the strength of a theory or a modelisation docs not really rely on the fact that the 
lithium observations are well reproduced. Other fits and signatures of the invoked processes 
have to be found to test the theory. 

We do not escape the rule : as pointed out in section 5.1, the diffusion coefficient 
parameters have been adjusted to obtain the best fit with the observations for lithium in 
the Hyades. However, once this has been done, there is no arbitrary parameter for the fit 
of the beryllium abundances, which is excellent. The results are thus very encouraging. 
Furthermore, we stress that the effective temperature for which the meridional circulation 
changes from one cell to two different cells is precisely that of the lithium dip, as already 
pointed out by CVZ92. This fact is independent from diffusion coefficient parameters. It 
can only slightly vary according to the prescription of the computations of the convective 
zone depth (here the mixing length theory, computed as for the Sun) . 

We claim that the dip feature is indeed related to this separation of the meridional 
circulation in two different cells. We have shown here that the two cells can be well dis- 
connected due to the process of "creeping paralysis" which occurs at the top of the deepest 
cell, thereby increasing the thickness of the frozen zone in-between. Consequences on the 
lithium abundances in sub-giants originating from stars hotter than the dip will be studied 
in a forthcoming paper (Theado, Vauclair and Dias do Nascimento 2003). 
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Fig. 5. — Lithium and Beryllium abundance variations at the surface of the stars for masses 
between 1.0 Mp,and 1.25 Up, 
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Fig. 6. — Lithium profiles inside the four cooler stars. The long dashed line represents the 
I.OM0 star, the dashed line : the I.IM© star, the solid line : the I.2M0 star and the dotted 
line : the 1.2M(^ star. 
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Fig. 7. — Beryllium profiles inside the four cooler stars with the same convention as the 
lithium profiles figure. 
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Fig. 8. — Evolution of /i-currents with radius inside a M = I.SMq star. The graphs show 
the variation with depth of both and E^^. 




Fig. 9. — Same figures as Fig. 8 for a M = l.SSM© star. 
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Fig. 11. — Same figures as Fig. 8 for a M = l.SM© star. 
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Fig. 12. — Meridional circulation in the stars with M>1.3OM0. In the less massive stars, 
the upper circulation loop disappears and the lower loop takes place in the whole radiative 
interior. 




Fig. 13. — Lithium and Beryllium abundance variations at the surface of the stars for masses 
between 1.3 Mp, and 1.5 M(7^ 
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Fig. 14. — Lithium profiles inside the four hotter stars. The long dashed-dotted line repre- 
sents the 1.3M0 star, the dashed line : the I.SSMq star, the solid line : the 1.4M0 star, the 
long dashed hne : the 1.45M0 star, the dotted hne : the I.SMq star and the long-dashed- 
dashed hne : the I.SSMq star. 
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Fig. 15. — Beryllium profiles inside the four hotter stars with the same convention as the 
lithium profiles figure. 
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Fig. 16. — Comparison between the computations of lithium (left) and beryllium (right) 
depletion and the observational values. The solid circles represents the models computed in 
these paper. Left figure : the open hexagons are observations from Boesgaard & Tripicco 
(1986), the open stars are from Boesgaard & Budge (1988), the open squares are from 
Thorburn ct al (1993) and the open circles arc observations from Boesgaard & Kind (2002), 
the inverted triangles denote upper limits. Right figure : the open circles are from Boesgaard 
& Kind (2002). 



